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Repulsive interaction helps superconductivity in fullerides 
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A repulsive interaction model of superconductivity (SC) is studied for tight-binding models with three- 
fold degenerate molecular orbitals. Taking a weak-coupling approach, we derive dimensionless coupling 
constants for various symmetries of SC pairs. In addition to anisotropic SC pairs, the s-wave pairing (A g ) 
can also be formed. With the purely repulsive interaction, however, the A g pair is not the most stable in 
both bcc and fee lattices. The most stable SC pair for the bec lattice has the T g symmetry, which is favored 
by a strongly nesting Fermi surface. In the fee lattice, various SC symmetries have comparable coupling 
strengths. With the electron-phonon interaction combined, it is likely that the A g pair becomes the most 
stable. 
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1. Introduction 

Alkali metal doped fullerides (A3C60) show not only super- 
conductivity with high transition temperature T c up to about 
40K, but also antiferromagnetism with A=Cs. It is known 
that CS3C60 can take two different cubic structures; A15 (bcc 
like) and fee lattices. 1 The conduction bands of A3C60 are 
composed of three-fold degenerate t\ u molecular orbitals of 
a Ceo- 2 Photoemission spectroscopy estimates the Coulomb 
repulsive interaction U in the range of 1-1.5 eV. 3 This magni- 
tude of U is substantially larger than the band width W (~ 0.5 
eV). 2 These estimates suggest that A3C60 is a strongly cor- 
related system. In fact CS3C60 is a Mott insulator at ambient 
pressure, 4 while other A3C60 with A = K, Rb are supercon- 
ductors at low temperature. 5-7 The insulating state of CS3C60 
changes to superconducting (SC) state under applied pressure, 
which is in line with other organic superconductor. 8 

It has been discussed that the SC state of A3C60 is driven 
mainly by the electron-phonon interaction. 9 The dynamic 
Jahn-Teller effect is invoked for reducing the strong Coulomb 
repulsion, and thereby stabilizing the ordinary s-wave state. 10 
In view of nearby presence of the antiferromagnetic state, 
however, the effective Coulomb repulsion should remain sig- 
nificant even in superconducting state. Then the basic ques- 
tion remains why the T c in fullerides is so high. 

The interplay of the Coulomb and electron-phonon interac- 
tions is a highly delicate matter in the strong-coupling region 
with high T c . As a possible step toward the complete under- 
standing, we take in this paper a complimentary approach, and 
study how the Coulomb repulsion helps superconductivity un- 
der the characteristic band structure with degenerate orbitals. 
We show that even with the repulsive interaction model, the 
s-wave pair can be formed. The situation is similar to that dis- 
cussed for iron pnictide compounds. 11 We further discuss (i) 
which symmetry is most stabilized in the presence of degen- 
erate orbitals, and (ii) how the difference between the Fermi 
surfaces of bcc and fee polymorphs influences the pairing. 

Among many energy bands in A3C60, we keep only the t\ u 
molecular orbitals located near the Fermi level. In order to 
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deal with the Coulomb repulsion, we follow and extend the 
weak-coupling theory for single orbital. 1215 Such extension 
has already been performed in the literature, especially for 
another multiband system Na A Co02 ■ yH^O. 16 Our treatment 
is similar to that in ref.16, but our motivation and conclusion 
is very different. 

In the following section, we present our model and discuss 
the perturbative approach in the presence of multi-orbitals. 
The explicit results for the coupling constants are presented 
in §3 for both bcc and fee lattices. Section 4 summarizes the 
result, and discusses possible relevance of our results to ac- 
tual fullerides. The technical details of calculation are given 
in Appendices A and B. 

2. Gap equation with multi-orbitals 

We focus on the three molecular orbitals of each with 
t\ u symmetry, and denote them as m = x, y, z. Although the ac- 
tual C60 molecules in A3C60 have two different orientations, 2 
we simply assume uniform orientation of C60 molecules. 
Then the space groups in our models become Bm3 (bcc) and 
Fm3 (fee). By going to the momentum space, we introduce 
the annihilation (creation) operator corner ( c l mcr ) f° r conduc- 
tion electrons with momentum k and spin cr. The band struc- 
ture is described by the effective tight-binding Hamiltonian 
as 



Wo = 2^ 2] e " m 

ka mn 



( k ^ C km<r Ck >">-> 



ka I 



Ckla = ^A ( "'\k)Ck lr 



(1) 



(2) 



(3) 



where / is the band index with the energy ei(k). The latter is 
obtained from diagonalization of the 3x3 matrix e mn {k) with 
use of Af^ik) as transformation from orbital (m) basis. De- 
tails of the hopping integrals involved in e mn {k) are explained 
in Appendix A. 

We consider the on-site Coulomb interactions such as intra- 
orbital Coulomb U, the inter-orbital Coulomb U', the Hund's 
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Fig. 1. Illustration of the KL interaction with m,n,m',n' being orbital in- 
dices. The diagrams in the upper part shows mean U and U%(k + k')U 
in the particle-hole channel. By exchanging the orbital labels in the lower 
line, we obtain KL interaction V mn j„> n >(k,k') in the particle-particle chan- 
nel as shown in the lower part. 



coupling J, and the pair-hopping J'. In deriving the pairing 
interaction originating from the Coulomb repulsion, we take 
the second-order perturbation theory according to ref.13, and 
refer to it as the Kohn-Luttinger (KL) interaction. 13 The KL 
interaction has been applied to several lattice models. 12,1415 
In the orbital degenerate system, the interaction depends on 
orbital indices m, n, m'n' ', and is given by 



V mn , m ;Ak, k') = [0+ U X (k + k')u\ 



(4) 



where x(k+ k') is the momentum-dependent static susceptibil- 
ity matrix. The Coulomb interaction U has matrix elements: 

Utnm,mm — Umn,mn ~ U 1 Umm,nn — J> Umn,nm — J 1 

n). The other components are set to zero. Figure 1 shows how 
to label the KL interaction with orbital indices. The interac- 
tion depends on the sum k + k' of the internal momenta of 
incident and scattered pairs. 

The matrix elements of £ are given by 



An 



,'n'(q) 



n2 h ' 1 

p,l,V 



(m) (p + q)A ( r'\p- 



f(e,(p + q)) - f(e,,{p)) 



(5) 



e t (p + q) - ev(p) 

where /(e) = l/(e /3(e_/j) + 1) is the Fermi distribution function 
with fi the chemical potential and /3 the inverse of temperature 
T. The bar in A ; (m) denotes the complex conjugate, and N is 
number of lattice points in the system. Note that^- sensitively 
reflects the shape of Fermi surface (FS), especially the nesting 
structure. Namely the KL interaction becomes strong if k + k' 
corresponds to a nesting vector, 

We derive the dimensionless SC coupling constant A in the 
following manner. Let us start with the linearized gap equa- 
tion given by 



N w 2ei<(k') 



(6) 



where A/(£) is the SC gap function with the band index /. Note 
that the interband pairing is not realized in the weak coupling 
limit. Accordingly the pairing interaction is also labeled by 
the band indices as 

Vw (k,k') = 

2 Af l \k)Af\-k)V mn , m ,,Ak, k')Af'\k')Af\-k'). (7) 

mn.m'n' 

We work with the weak coupling limit where A/(£) is siz- 
able only near the FS. Since T c is much smaller than the char- 
acteristic energy of we can neglect the ^-dependence of the 
effective interaction V mn ,m'n'(k, k'), and put T = in eq.(5). 
Then, we can rewrite eq.(6) as 



A,(k) - -i V Viv{k,k')8(e k , v )A v {k' 



r*Lt) c 
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de 



tanh(e/2r c ) 



(8) 



where a> c is the cut-off energy. By changing the summation 
over k by the surface integral over the FS, we obtain the eigen- 
value equation for A as 

where dS m denotes a FS element of the band / with v^(£, I) 
being its Fermi velocity, and 

tanh(e/2T c ) 
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Equation (9) shows that A is independent of co c , while eq.(10) 
shows that T c ~ <d c e\p(-l/A) does depend on ai c . In the fol- 
lowing, we concentrate on A rather than T c . We evaluate the 
surface integral by using the tetrahedron method, 17 and the 
eigenvalue equation (9) is solved by the power method. In our 
calculation, more than 10 4 elements on FS are used. 

In general, a SC order parameter involves all conduction 
bands with a component A/(£) in each band /. The largest 
eigenvalue leads to the highest T c and its eigenvector de- 
termines the symmetry of the SC state. 18 Since each Ceo 
molecule does not have the four-fold rotation axis, the SC 
state belongs to the irreducible representation of point group 
7 1 /,. 19 Then, some gap functions in the cubic (9;, group with 
higher symmetry are mixed under the T/ t symmetry. 20 For 
the spin singlet Cooper pair, there are three representations: 
A g , E g , T g . Imposing these symmetries on A/(£), we reduce 
the £-space integration in eq (9) to a smaller part of the BZ. 
In this way we obtain eigenvalue equations for each symme- 
try. More details of the reduction procedure is explained in 
Appendix B. In order to avoid the O(U) term in the KL inter- 
action, a nodal SC state is favored in a single band model. 12 
However, since A/(£) can change the phase factor at each 
band, even a fully gapped SC state can avoid the O(U) term 
in the multiband model. Note that similar situation has been 
discussed in iron pnictide compounds. 11 



J. Phys. Soc. Jpn. 



Full Paper 



Author Name 3 



Band2 



Table I. The SC coupling constant A for each irreducible presentation (IR). 
The unit of U is eV, which means that /! takes the numerical value in the ta- 
ble with U = leV, although the relevant magnitude for U is much smaller. 



IR 


Degeneracy 


A/U 2 (bcc) 


A/U 2 (fee) 


A s 


1 


0.2358 


0.2016 


E g 


2 


0.1349 


0.2926 


T g 


3 


0.6287 


0.2552 


T u 


3 


0.2860 


0.1869 
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Fig. 2. The conduction bands derived in the tight-binding bcc model. 

3. Numerical Results 

3.1 Choice of parameters 

We have derived the maximum A for each symmetry of 
singlet pairing with fee and bcc crystal structures. Since the 
triplet pair is less interesting for fullerides, we have only de- 
rived A with T u symmetry. In the following, we consider only 
the case of three conduction electrons per Cgo- We fix the band 
width W to 0.5 eV. Concerning the exchange interaction, it 
has been argued that the Jahn-Teller interaction can change 
the sign of the Hund's rule coupling. 21 In view of such uncer- 
tainty, we simply put J — J' — 0, and U — U' . Then we obtain 
the simple scaling A cc U 2 , provided only the 0(U 2 ) term in 
eq.(4) is active. The values of A computed for each symme- 
try are summarized in Table I. Note that the theory remains 
accurate only for U much smaller than the band width W. 

3.2 bcc model 

We have fitted the tight-binding parameters so as to repro- 
duce the energy bands derived by a more elaborate method. 22 
With only the nearest- and next-nearest hoppings included, 
a satisfactory fitting can be performed, as explained in Ap- 
pendix A. Figure 2 shows the calculated result for conduction 
bands in the bcc model. The resultant FS is shown in Fig. 3. 
The FS in band 1 is shaped like a cube, while the FS's in band 
2 and band 3 are not closed the first BZ. As a result the FS 
consists of six sheets that are very flat in band 3, and less flat 
in band 2. If we combine six sheets at each band beyond the 
first BZ, the FS in band 2 and band 3 also looks like a cube 
centered on H point in the BZ. Note that a perfect cube of 
the FS is obtained in a half-filled single band model with only 
the nearest neighbor hopping in the bcc lattice. This feature is 
weakened by band mixing, but remains to some extent in the 
present bcc model. 

The strong nesting in the FS appears in the momentum de- 
pendence of the static susceptibility xiq)- Figure 4 shows the 
main components Xmn,mn- The other components Xmm,nn an d 



Fig. 3. The FS in the the bcc model. The center of the BZ is T. See text 
for comment on the shape. The FS in band 1 is shaped like a cube. FSs in 
band 2 and band 3 have six sheets which are almost flat. If we combine 6 
sheets at each bands, FS in band 2 and band 3 also look like cube which is 
centered in H point. In single band bcc model, FS becomes a cubic shape 
perfectly. 
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Fig. 4. Momentum dependence of the main components Xmn.mn(<l) (»!, n = 
x, y, z). Location of high symmetry q points is illustrated in Fig. 3. The 
enhancement near H point is related to antiferromagnetic fluctuation in the 
bcc model. Note that summation of Xmm,m °ver m and n recovers the Tj, 
symmetry. 



Xmn,nm with m + n are smaller by an order of magnitude 
than Xmn,mn- Those components with three different indices 
such as Xxy,zx are zero since each molecular orbital is odd un- 
der space inversion. In Fig. 4, the enhanced response around 
q = Q H = (2tt, 0,0), (0,2^,0), (0,0, 2tt) is related to insta- 
bility toward antiferromagnetic (AF) order. This instability is 
strong in the bcc lattice with bipartite structure. In the single 
band model with perfect nesting at half filling, the AF suscep- 
tibility x(Qh) diverges at Q H . 

Using the susceptibility thus derived, we calculate the SC 
coupling constant. The largest eigenvalue of A in eq.(9) for 
each symmetry are shown in Table I. In the bcc model, the 
T g symmetry of the pair is most favored. The eigenstate has 
the three-fold degeneracy with components {d xy , d yz , d zx ). In 
order to understand the stability of the T g pairing, we show in 
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Fig. 5. Results for the gap function A/(fc) with d xy symmetry projected onto 
the plane k z = in the bcc model. The order parameter with T g symme- 
try consists of three components d xy ,d yz ,d zx . The dashed lines show the 
boundary of bcc BZ, and the dotted line shows the completely nesting FS 
of the single band model with only the nearest-neighbor hopping. The d xv 
state has nodal points at k x = and k y = 0, which extends to lines by 
including finite k z . 



Fig. 5 the gap function A/(k) of the d xy state along the FS in 
the k z = plane in BZ. The numerical values shown by color 
mean the amplitude with normalization over the whole Fermi 
surface. The red and blue parts represent positive and negative 
A/(£), respectively, that can be taken real. The boundary be- 
tween red part and blue part becomes a node. For comparison 
the completely nested FS in the single-band case is shown by 
dotted lines. The band 2 and band 3 satisfy the nesting condi- 
tion reasonably well. 

It is known in the single band bcc model 12 that the SC state 
with T g symmetry is most stabilized. To understand this as 
well as our result for the multiple bands, suppose in Fig. 5 
that a SC pair with red region of k is scattered to k' in the 
blue region in the same band. With the d xy symmetry, there 
is substantial combination to satisfy k + k' ~ Qh- Hence, the 
effective interaction Vu>(k, k') becomes effective for c/ lv state. 
Similar situation occurs to other members d yz , d zx of the T g 
symmetry. Note that the gap functions of all bands have the 
same sign for each quadrant of BZ in Fig. 5. This is because 
each band has comparable size of FS, and compose almost 
degenerate bands. 

We next move to the pairing with the A g symmetry which 
is less stable in the present model, but which seems most rel- 
evant to actual fullerides. Figure 6 shows the gap function 
A/(k) for each band. Different signs of the (real) gap func- 
tions work to cancel the O(U) part of the KL interaction. We 
emphasize such cancellation requires nodes along the FS in 
the single band model. In our case, the band 2 has nodes as 
in the single-band model, while band 1 and 3 cancel the 0{U) 
without nodes, and hence having less cost in the kinetic en- 
ergy. 



Fig. 6. Illustration of the gap function A/(fc) with A g symmetry projected 
onto the plane k z = in the bcc model. Note that the gap functions of 
bands 1 and 3 have opposite signs, while the band 2 has nodes keeping the 
A g symmetry. 




Fig. 7. The conduction bands in the tight-binding fee model. 



3.3 fee model 

The tight-binding fit of the band structure of CS3C60 is 
available in literature. 2 We have used the same fitting parame- 
ters and reproduced the previous result. For completeness, we 
show our result for the fee model in Fig. 7. Figure 8 shows the 
FS in the fee model. Although the conduction bands consist 
of three branches, there are only two FS's. It is clear in band 
2 that there is no four-fold symmetry in the FS, which is how- 
ever consistent with the local 7 1 /, symmetry of Ceo- Note the 
presence of the three-fold symmetry around the [111] axis. 

Figure 9 shows the momentum dependence of the main 
components Xmn,mn- The other components are negligible as 
in the bcc model. The average of Xmm,mm is shown by right 
blue line, which does not have a particular wave vector for 
enhanced behavior. This corresponds to a difficulty to realize 
AF order in fee lattice which has a geometrical frustration. 
Using this susceptibility, we calculate the SC coupling con- 
stant. In Table I, the SC state with E g symmetry is the most 
stable in the fee model, but the difference with other symme- 
try is not large. 

Figure 10 shows the gap function Ai(k) with the E g sym- 
metry in the fee model in the k, = plane of BZ. In con- 
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Fig. 8. Two sheets of the FS in the fee model. High symmetry points in the 
first BZ are shown. See text for comments on the shapes. 




Fig. 9. The momentum dependence of the main components Xmn,mn (m,n = 
x, y, z) together with the average. High symmetry points are illustrated in 
Fig. 8. From X point to W point, some ^„„„„„„ becomes large, but the aver- 
age of Xmm ,mm has no special peak. 




k x /ji k x /7i 



Fig. 10. Illustration of the gap function A;(t) with E g symmetry in the fee 
model. The dotted lines show the boundary of fee BZ. The nodes in band 
2 has no four-fold symmetry because of the T/, symmetry of each Cgo 
molecule. 

trast with the (9/, group, which has d x 2_ Y i and dj tZ i_,2 states as 
eigenfunctions, the eigenfunctions in the T/, symmetry have 
no four-fold symmetry. In our result, A^k) corresponding to 
band 1 is close to d x %-^ state, while ^(k) is rather different 
from the E g state in the 0/ t group. The node structure of Ai(k) 
in Eg symmetry is complicated. In band 1, the nodal structure 
is rotated by tt/4- from the corresponding band in the bec case. 
On the other hand, in band 2, the node appears not only along 
the FS but also in disconnected parts along k y = +2n of the 
FS. 

Finally we present the results for the A g pairing in the fee 
model. The coupling constant A as given in Table 1 is almost 
the same as in the bec model. Figure 1 1 shows the gap func- 
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Fig. 1 1 . Illustration of the gap function A/(k) with A g symmetry in the k. = 
plane in the fee model. The two sheets of FS have opposite signs of Aj(ifc). 

tion in the k z — plane of the BZ. There is no node for each 
FS in accordance with the s-wave like character. However, the 
sign of A[(k) is opposite between bands 1 and 2. In this way, 
the first-order Coulomb repulsion is canceled. 

4. Summary and Discussion 

We have studied the repulsive interaction model for super- 
conductivity in A3C60 By using the second order perturbation 
theory with respect to the Coulomb repulsion, we calculate 
the SC coupling constant A for various SC symmetries and 
for both bec and fee lattices. By the nesting property of the 
FS in the bec model, the T g symmetry is the most favorable 
within the repulsive interaction model. We have further shown 
that the stable pair with the fully symmetric A g can also be 
formed by purely repulsive interaction. This stability essen- 
tially requires the presence of multiple conduction bands. 

Let us now discuss possible relevance of our results to ac- 
tual superconductivity in fullerenes, especially in CS3C60 with 
the highest T c , The most fundamental question is why the T c 
is so high. Since the system is close to the Mott transition, one 
might naively guess that the SC is caused by a non-phonon 
mechanism as in cuprate superconductors. However, various 
experimental evidences point to the fully symmetric s-wave 
SC being realized. In the conventional theory, the s-wave SC 
is unfavorable in the presence of strong Coulomb repulsion. 

According to our results in this paper, we propose the fol- 
lowing scenario for the SC in CS3C60: What is responsible for 
the high T c is the cooperation, rather than competition, be- 
tween the Coulomb repulsion and the Jahn-Teller phonons. 
The A g pair will be most favorable for such cooperation, 
which is fully gapped in the fee model. In the bec model, on 
the other hand, one of the three bands has four-fold nodes in 
the gap function. Hence our repulsive interaction model pre- 
dicts gapless s-wave superconductivity in the bec system. It is 
not clear whether the nodes remain in the presence of Jahn- 
Teller phonons. Theoretically, simultaneous account of Jahn- 
Teller phonons together with Coulomb repulsion requires a 
new scheme which is not available at present. The main 
difficulty is that the KL interaction is highly non-local and 
hence is beyond the scope of the dynamical mean-field theory 
(DMFT). We have of course recognized that the DMFT can 
powerfully deal with local Coulomb correlation together with 
Jahn-Teller phonons. 10 Various available schemes in the mo- 
mentum space, on the other hand, are not reliable enough to 
deal with the case of strong Coulomb repulsion. Namely, SC 



6 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



in fullerenes provides a challenging ground to construct a new 
theoretical scheme. We hope to contribute to further develop- 
ment of the theory in the near future, and test the scenario 
mentioned above. 
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Appendix A: Hopping integrals for t lu molecule orbitals 

We explain how the parameters are chosen in the effec- 
tive tight binding models for A3C60. Using the molecular or- 
bitals specified by m — x, y, z, we consider the hopping inte- 
grals t„ m {a,l3,y) between Ceo molecules with the directional 
cosines a,/3,y. Then the effective tight binding Hamiltonian 
is written as 



rameters for t n along (110): 



Wo=^]^ t mn (a,p, y) {c] ma .c jmT + h.c) 

<i,p 

-z 



(A-l) 



o~,mn 


Ckxcr 


t / 








V 



e z ,(k) 



e zv (k) 



(A-2) 



Ckx 
Ckzo-j 

is a summation over the nearest and the next 



e vz (k) 

€zz(k)) 



where £<,-,> 

nearest neighbor sites. The parameters t m „(a,/3, y) are fitted to 
be consistent with first principle band calculations. 2 ' 22 

A. 1 bcc model 

We keep only the the nearest-neighbor hopping i„ and the 
next-nearest-neighbor one f„„. the bcc model. Among eight 
nearest neighbors and six next-nearest neighbors, we take as 
representative members, a nearest neighbor along (1,1,1) and 
a next-nearest neighbor along (1,0, 0). The hopping integrals 
are parameterized as 



<X 
Y 



Y 
X 
Y 



Y\ 

Y 

X 



t tin — 



(Z 







w 








w 



where three orbitals of 1 1„ have been taken as the basis of the 
matrix. The hopping integrals to other equivalent neighbors 
are determined by symmetry. Fitting to a first-principle calcu- 
lation, 22 we choose the following values: 



X = 0.026, 
Z = 0.016, 



Y = -0.015, 
W = -0.008, 



W = -0.025. (A-3) 



The inequality W + W comes from the absence of the four- 
fold symmetry in T n , which is beyond the Slater-Koster pa- 
rameterization for p orbitals. 23 

Note also that actual CS3C60 takes the A15 structure where 
Ceo molecules have different orientations in the body-center 
and corner positions. Our form of i„ assumes, following 
ref.22, a simplified form that C60 molecules are all equivalent. 

A.2 fee model 

The fitting parameters for fee model have been introduced 
already by Gelfand et al. 2 Following their results, we keep 
only the nearest-neighbor hopping and use the following pa- 





X 


Y 


0" 


in = 


Y 


X' 







lo 





z) 



with 

X = 0.0083, X' = 0.0336, Y = -0.0198, Z = -0.0191. 

The inequality X + X' reflects the absence of four-fold sym- 
metry in each C60 molecule. The actual C60 molecules in the 
fee A3C60 take one of two orientations in a random way. 24 We 
assume for simplicity the same Ceo orientations as in the bcc 
model. 

Appendix B: Use of Symmetry in Deriving Eigenvalues 

We describe how to utilize the symmetry of A/(k) in deriv- 
ing the eigenvalue. For an irreducible representation a, eq.(9) 
takes the form: 

1 v r 



A a Ai a (k) 



(2n? Z-jJ eil , =eF v F (k',l') 



w ^-r£,(*,*')A,, a (*'), 



(B-l) 



where F^,(£, k') is the effective interaction. It is possible to 
reduce the range of integration to k' x , k' v , k[ > by using the 
symmetry operation. Table B-l shows the relation between 
each symmetry operator O and irreducible representations of 
7),. The basic symmetry operators change the wave number k 
and the phase factor of the gap function. 

Note that we do not necessarily use all operators for reduc- 
ing the BZ integration. For example, since C3 and C 2 change 
d„ basis function to other basis function, they are not utilized 
in the case of T g . In this paper, we choose A g , Eg, and d xy as 
basis functions. Then, T" t ,{k, k') for each symmetry are given 
in terms of the KL interaction Vw(k, k') by 

rjf(t, k') = [l+ o-Ak')] ■ [1 + o- y (k')} • [1 + cr z (k')] 

x [1 +C 3 (k') + Cj(k'Wii'(k,k'), (B-2) 

T%(k, *') = [! + o-Ak')] • [1 + or y (k')] • [1 + cr z (k')] 

x [1 + coC 3 (k') + to 2 Cl(k')]Vu>(k,k'), (B-3) 

r£(* t k') = [l- o-Ak')] ■ [1 + o-Ak')] • [1 + o-Ak')] 

xV„>(k,k'), (B-4) 

with a) = exp(27r//3). 

We obtain the maximum eigenvalue A a by using eq. (B-l) 
for each basis function. For a degenerate representation such 
as T g , three functions such as d xy ,d yz and d zx may enter into 
the process of iteration in the power method. However, we 
actually need only d„-type functions as shown below. 

Let us consider a matrix element 



^) 6 ?1Y' 



dS 1 



dS 



k'V 



{2nf^J J v F (k,Qv F (k',l') 

x (A%\kjfv u >(k,k')A%(k'), (B-5) 
where m, n represent either xy, yz or zx in T g . In our multiband 
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Table Bl. Eigenvalues of the basic symmetry operation O for each irre- 
ducible representation of the gap function A;. Degenerate representations 
such as E g and T g can have diagonal basis set by proper choice given in 
the Table. Namely, d„,, d yz , and d- x are the basis functions of T g symmetry. 
Each operator O also changes the wave vector k as shown in the second 
column, and the phase factor as shown from the third columns. For ex- 
ample, the operator <x v changes the gap function with d„ : symmetry as: 

Al(-k X ,ky,k ! ) = (-l)XAl(k X ,ky,k Z ), 






Ok 




4 


c « 


dxy 


dy 7 


d z . 


E 


k,K-> ky, Kz 


i 


i 


1 


1 


1 


1 


(fx 


k>xt ky, k-7 


i 


i 


1 


-1 


1 


-1 


(Ty 




i 


i 


1 


-1 


-1 


1 


°~: 


ky, ~k z 


i 


i 


1 


1 


-1 


-1 




ky • k.-'y * k ^ 


i 


OJ 


oj- 








c 2 


k-7 . k^, ky 


i 


0) 


(Jj 









model, AT'(k) and Vw(k, k') satisfy the relations: 

a x (k) ■ A%\k) = A^\-k x ,k y ,k z ) = s m ■ A^(k x ,k v ,k z ), 

(B-6) 

a x {k) ■ a x {k') ■ V u .(k, k') = V,K-k x , k y , k z , -k' x , k[„ k' z ) 

= V u ,(k x ,ky,k z ,k' x ,k' y ,k' z ), (B-7) 

where s m is a phase factor of o~ x . Now we restrict the region 
of integration to k x > 0, k\ > 0. Then, /„,„ is rewritten as 

1 yi ff dS k i dS k 'i> 
"'" ~ (2;r) 6 Zj J J VF {k,l)v F (k',l') 

x {A^(k))*V w {k,k')Af a {k'), (B-8) 

V u ,(k, k') = [(l + + (s m + s n )o- x (k)]Vu'{k,k'). 

(B-9) 

When we put m = xy and n = yz, for example, we obtain 
s,„ = -1, s n = 1. Then, Vw{k,k') and hence /,„„ become 
zero. Also for other cases including E g , we generally find that 
htm = for m + n. Therefore we need only a single basis 
function in to obtain the maximum eigenvalue A a for each a. 
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